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ABSTRACT 


Laser and camera data taken during the International Satellite 
Geodesy Experiment (ISAGEX) have been used in dynamical 
solutions to obtain center-of-mass coordinates for the Astro- 
Soviet camera sites at Helwan, Egypt and Oulan Bator, Mon- 
golia and East European camera sites at Potsdam, German 
Democratic Republic and Ondrejov, Czechoslovakia. The 
results are accurate to about 20m in each coordinate. The 
orbit of PEOLE (i = 15°) has also been determined from 
ISAGEX data and mean Kepler elements suitable for geo- 
dynamic investigations are presented. 


“"Presently with the National Oceanic and Atmospheric Administration/National Ocean Survey, 
Rockville, Maryland. 
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GEODETIC RESULTS FROM 
ISAGEX DATA 


INTRODUCTION 

'Hie ISAGEX experiment was initiated by the French Centre National d'Etudes 
Spatiales (CNES) in the autumn of 1969 through a proposal foF anlnte rnational 
laser and photographic campaign on satellites equipped with laser reflectors. 

The ISAGEX experiment was endorsed by the COSPAR XHIth General Assembly, 
Leningrad, 1970, The objective of the program was to collect a set of homoge- 
neous and well distributed precise laser and camera satellite observations for 
the purposes of dynamic and geometric geodesy. The experiment involved 
seventeen countries and over fifty tracking stations. The data gathering por- 
tion of the experiment extended from December 1970 to September 1971. 

The results of a dynamical solution which combined ISAGEX laser data with the 
optical and laser tracking data recorded during the National Geodetic Satellite 
Program (NGSP) and the CNES/SAO 1968 Observing Program have been re- 
ported previously by Marsh, et al. (1973). This solution provided coordinates 
for over 70 globally distributed tracking stations. Also very preliminary analy- 
ses of Astro-Soviet and East European camera data received as of May 1973 
were reported by Marsh, Douglas and Klosko (1973b). The present paper repre- 
sents a continuation of the analyses of the ISAGEX data using the laser and 
optical data available as of May 1974 to refine coordinates for the camera 
stations. 


DYNAMICAL SOLUTIONS FOR STATION COORDINATES 

The orbital and station coordinate recovery solutions were derived through the 
use of the GEODYN program (Martin, 1972) on the GSFC IBM 360/95 computer. 
GEODYN is a multiple arc, multiple satellite orbit and geodetic parameter 
estimation system based upon Cowell-type numerical integration techniques. 
Model parameters included luni-solar gravitational perturbations, solar radia- 
tion pressure, the Jacchia Model Atmosphere (1965, 1968, 1970, 1971) for drag 
computation, BIH polar motion and UTI data, and the GEM-1 gravity model 
(Lerch, et al. 1972) modified by the use of the resonant coefficients from the 
SAO 1969 S.E, II (Gaposchkin and Lambeck, 1970). This gravity model was used 
to insure consistency with the GSFC 1973 solution. 

In the present effort, we have been primarily concerned with computation of 
orbital residuals from the Astro-Soviet and Eastern European data and a 
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dynamical adjustment of the coordinates for these sites since data from these 
stations have not previously been available. The stations are: 


ASTRO-SOVIET EAST EUROPEAN 

Oulan Bator, Mongolia Ondrejov, Czechoslovakia 

Helwan, Egypt Potsdam, German Democratic Republic 

The Astro-Soviet instruments were the AFU-75 tracking cameras and the East 
European stations used the Zeiss SBG camera. Figure 1 shows the locations of 
these stations, along with the other stations which participated in the ISAGEX. 
Table 1 presents a summary of the data used from these four stations, along 
with the RMS values of the residuals in the solutions. Data was contributed by 
several other Astro-Soviet and East European camera stations but were not 
analyzed since sufficient numbers of observations were not available for dy- 
namical orbit computations. However, the data could be used in geometric 
analyses. 

These data were preprocessed with corrections applied in accordance with the 
information provided in ISAGEX Report No. 16, "Data Handling Booklet" 

(Brachet 1973). In the case of the Astro-Soviet data, we applied corrections for 
annual and diurnal aberration and parallactic refraction; precession and nutation 
were applied to convert from the reference system with a mean equator and 


Table 1 

S umm ary of Astro- Soviet and East Europe Camera Data 
Used in Dynamical Computations 


Station 

Right Ascension* (a) 

Declination ( S) 



No. 

R.M.S. 

No. 

R.M.S. 

Oulan Bator, Mongolia 

27 

3.0 arc sec 

29 

4,2 arc sec 

Helwan, Egypt 

46 

2,9 arc sec 

47 

2.8 arc sec 

Ondrejov, Czechoslovakia 

32 

3.5 arc sec 

32 

6.0 arc sec 

Potsdam, German 
Democratic Republic 

25 

4.0 arc sec 

26 

4.6 arc sec 


* Right Ascension Residuals have been multiplied by cos 8. 
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GL GRGS LASER W WEST EUROPEAN CAMERA 


Figure 1. Locations of Stations Which Collected Data During the International Satellite Geodesy Experiment 







equinox of 1950.0 to a true of date system. The passive observation time tags 
were converted from USSR UT1 time to U.S. Naval Observatory (USNO) UTC 
time. 

Corrections for annual and diurnal aberration were already applied to the data 
from Potsdam (1181) and Ondrejov (1147). We applied parallactic refraction to 
the Ondrejov data. The passive observation times for these data were converted 
from_the BIH UTC system to the USNO UTC system. 

In this solution, coordinates for the ISAGEX lasers and other cameras were held 
fixed at the values derived earlier by Marsh, et al., (1973). Thus, these new 
coordinates were determined from combined laser/ optical reference orbits. 

In the interest of efficient use of the computer, and since little overlap existed 
in the data, separate multi-arc solutions were made for each of the 4 site loca- 
tions. Table 2 presents the arcs used in these solutions. Note that the solution 
for Helwan consisted of GEOS-1 data alone. GEOS-2 data could not be used be- 
cause of large inconsistencies. If this problem could be solved, the solution 
would probably be significantly improved. In contrast, only GEOS-2 data 
(flashes) could be used for the Ondrejov solution, since there was a problem of 
timing system fluctuations as large as 7ms for passive satellites at this station. 


EVALUATION OF STATION COORDINATES 

Table 3 presents the results obtained from this investigation and those provi- 
sional values given in ISAGEX Report No. 7 (Brachet, 1970). Differences be- 
tween the provisionally available coordinates and the new values range from a 
few tens of meters to a few hundred meters. 

The largest error source in this analysis is due to uncertainty in the modeling 
of the Earth's gravity field. As a means of assessing the magnitude of this 
error, the solutions were computed using the SAO 1969 (Gaposchkin and Lam- 
beck, 1970) gravity model as well as the GSFC GEM-1 model mentioned earlier. 
The differences in the adjusted latitude and longitude coordinates were less than 
20 meters in every case, while the height coordinate differences were less than 
10 meters. Since the GEM-1 model gave a lower (by almost a factor of two) 
residual rms it is concluded that the GEM-1 solution is more accurate. Thus, 
the coordinate errors due to gravity model uncertainty are probably less than 
the differences in the two sets of coordinates and are assessed as less than 
20 meters in latitude and longitude and less than 10 meters in height. 

An initial height comparison at Oulan Bator indicated a difference on the order 
of 1400 meters between the dynamically derived value and the value of 175 m 
provided in the ISAGEX Report No. 7 (Brachet, 1970). Further checking with the 
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Table 2 

Orbital Arcs Used in. Dynamic Solutions 



Satellite 

Begin 

End 





hr. 

“min. 



yr. 

mo. 

day 

nr. 

min. 

yr; 

mo, 

day 

Oulan Bator 

Arc 1 

GEOS 1 

71 

02 

23 

02 

15 

71 

02 

25 

22 

14 

Arc 2 

GEOS 1 

71 

02 

26 

02 

31 

71 

02 

28 

22 

28 

Arc 3 

GEOS 2 

71 

04 

13 

18 

40 

71 

04 

14 

21 

47 

Arc 4 

GEOS 2 

71 

04 

16 

07 

28 

71 

04 

19 

04 

12 

Potsdam 

Arc 1 

GEOS 1 

71 

03 

02 

00 

40 

71 

03 

06 

20 

54 

Arc 2 

GEOS 2 

71 

04 

03 

07 

00 

71 

04 

05 

05 

09 

Arc 3 

GEOS 1 

71 

07 

01 

02 

29 

71 

07 

03 

02. 

53 

Ondrejov 

Arc 1 

GEOS 2 

71 

04 

01 

01 

48 

71 

04 

04 

19 

31 

Arc 2 

GEOS 2 

71 

04 

13 

02 

14 

71 

04 

15 

23 

22 

Arc 3 

GEOS 2 

71 

04 

19 

04 

05 

71 

04 

21 

08 

35 

Arc 4 

GEOS 2 

71 

04 

29 

01 

51 

71 

05 

02 

20 

32 

Arc 5 

GEOS 2 

71 

08 

13 

01 

11 

71 

08 

16 

00 

40 

| Helwan 

Arc 1 

GEOS 1 

71 

07 

05 

00 

20 

71 

07 

08 

01 

09 

Arc 2 

GEOS 1 

71 

07 

11 

22 

44 

71 

07 

14 

03 

00 

Arc 3 

GEOS 1 

71 

07 

20 

21 

41 

71 

07 

22 

22 

02 

Arc 4 

GEOS 1 

71 

07 

23 

19 

31 

71 

07 

25 

23 

55 

Arc 5 

GEOS 1 

71 

07 

26 

20 

01 

71 

07 

30 

20 

32 
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Table 3 


Astro-Soviet and East European Station Coordinate Values 


Station 


Latitude 

Longitude (E.) 

Height 


(Deg, 

Min, Sec) 

(Deg, Min, Sec) 

(Meters) 

Oulan Bator, 

New Value 

47 

m 

55.1 

107 

3 

10.7 

1 

Mongolia 

Old Value* 

47 

E9 

56 

107 

3 

00 



Difference 



.9 



10.7 

35 

Helwan, Egypt 

New Value 

29 


31,6 

31 

20 

32.5 

97 


Old Value* 

29 

51 

31.1 

31 

20 

28.05 

120 


Difference 



.5 



4.45 

23 

Ondrejov, 

New Value 

49 

55 

15.8 

14 

47 

55.4 

525 

Czechoslovakia 

Old Value* 

49 

55 

19.4 

14 

48 

3.9 

537 

1 

Difference 



3.6 



8.5 

12 

Potsdam, GDR 

New Value 

52 

22 

49.7 

13 

3 

54.1 

147 


Old Value* 

52 

22 

51.4 

13 

3 

58.8 

136 


Difference 



1.7 



4.7 

11 

Coordinates are referred to an ellipsoid with: 




a e = 6378155 meters, l/f = 298.255 




*From Brachet, 1970. 


Defense Mapping Agency/ Topographic Command (Niska, 1974) revealed that the 
MSL height should be 1550 m. Thus using a geoid height of -45 meters (Marsh 
and Vincent, 1974) results in an ellipsoidal height of 1505 m or a difference of 
35 meters with respect to the dynamic value. 

An additional comparison was possible for the height coordinate at Potsdam 
since a value for the surveyed height above mean sea level (MSL) was provided. 
To this value of 109 meters, the geoid height of 27 meters from the GEM-6 de- 
tailed gravimetric geoid (Marsh and Vincent, 1974) was added resulting in an 
implied height above the reference ellipsoid of 136 meters, 11 meters less than 
the value recovered dynamically. This level of agreement is considered good 
since the 11 meter difference includes errors in the dynamical recovery, MSL 
survey errors, and errors in the geoid. 

The quality of simultaneous data is important to consider. Figures 2 and 3 show 
the residuals in right ascension times cosine declination (a cos 8) and declina- 
tion (S) for typical simultaneous passes of Helwan with Dionysos, Greece and 
with Debra Zeit (DEZEIT), Ethiopia. The consistency is generally very good. 
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JULY 6, 1971 

22 h 31 m UTC 
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Figure 2. 

Residuals from Simultaneous Observation of GEOS-1 by Cameras 
at Dionysos, Greece, {G), and Helwan, Egypt, (H) 


A a COS 8 (ARC SECS) 


00 



SECONDS FROM 
JULY 28, 1971 

17 h 56 m UTC 



* EDITED FROM SOLUTION 


Figure 3. Residuals from Simultaneous Observation of GEOS-1 by Cameras 
at Helwan, Egypt, (H) and Dezeit, Ethiopia, (D) 




except for the a cos S residuals in Figure 3. The inconsistency is about 3 arc 
seconds between the centroid of the Helwan residuals and the centroid of the 
DEZEIT residuals. 

Figures 4 and 5 show residuals for simultaneous tracking of Ondrejov with Haute 
Provence and/or Zimmerwald. The consistency is excellent, small average 
differences of only the order of one arc second appearing. 


MEAN ELEMENTS OF PEOLE* 

During 1971, osculating orbital elements for PEOLE were determined at 4-6 day 
intervals to provide prediction data for the GSFC laser at Guam as part of the 
International Satellite Geodesy Experiment. The orbits were determined almost 
entirely from CNES and NASA Minitrack data, with occasional passes of GSFC 
laser data from Guam, CNES laser data from Dakar, Senegal and SAO laser data 
from Arequipa, Peru and Natal, Brazil also employed. Because of the low in- 
clination of the orbit and resultant diminished geopotential perturbations, the 
orbits determined from the Minitrack data alone are surprisingly precise. The 
uncertainty of the elements is much lower than the amplitudes of the orbit vari- 
ations of interest. 

The orbit of PEOLE is especially interesting for atmospheric and geopotential 
investigations. The mean decay rate of the semi-major axis is about 30m/day, 
but there is considerable variation over the year. The perturbations due to odd 
zonal harmonics are also very large. Solar tidal perturbations may also be 
perceptible in the inclination. 

We have prepared mean orbital elements of PEOLE from the osculating elements 
by the combined analytic-numerical technique developed by Douglas, Marsh, and 
Mullins (1972). The method involves the analytical removal of all large short 
periodic perturbations followed by a numerical averaging over 1 day to remove 
remaining small high frequency effects. Thus all long periodic and secular 
effects remain in the mean elements. 

Figure 6 shows the mean semi-major axes obtained from the orbital fits. The 
SAO 1969 Standard Earth (Gaposchkin and Lambeck, 1970) gravity model was 
used in the orbit determinations. To get the most accurate orbits possible, the 
drag coefficient was adjusted in each orbital arc. When the mean elements were 
computed, the appropriate drag coefficient was used for each computation. Thus 
the mean elements given here are an accurate portrayal of the evolution of the 
orbit. 


^Presented to the Symposium on Satellite and Terrestial Triangulation, Graz, Austria, May 1972. 
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Figure 4. Residuals from Simultaneous Observation* of GEOS-2 by Cameras at Haute Provence, France, (P) 

Zimmerwald, Switzerland, (Z), and Ondrejov, Czechoslovakia, (O) 
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Figure 5. Residuals from Simultaneous Observation* of GEOS-2 by Cameras at Ondrejov 
Czechoslovakia, (O) and Zimmervvald, Switzerland, (Z) 
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Figure 6. Mean Semi-Major Axis PEOLE-ISAGEX 


As is evident in Figure 6, the PEOLE orbit decays rapidly, up to 1 Km per 
month. This can be a serious matter for geodetic investigations unless the drag 
perturbation can be accurately modeled. Preliminary results show that if the 
drag coefficient C D is allowed to adjust and the atmospheric density is modeled 
by the Jacchia model atmosphere (Jacchia 1965, 1968), the unmodeled part of 
the decay is only a few hundred meters out of a total decay of about 12 Km 
during the year. 

Figure 7 shows the mean eccentricities of the PEOLE orbit. Clearly visible 
are the secular decay due to atmospheric drag and the approximately 28 day 
oscillation due to the odd zonal harmonics. Figure 8 shows the residual eccen- 
tricity for the first six months of 1971 after removal of atmospheric drag effects 
by the Jacchia model, zonal harmonic effects by the SAO 1969 Standard Earth, 
and luni-solar perturbations. About 80% of the odd zonal harmonic effects have 
been accounted for, but a discernable unmodeled 28 day effect remains. The 
precision of the mean eccentricities is probably better than one part in 10 s , 
which is consistent with the observed dispersion in the semi-major axis. 

Table 4 presents all of the mean elements of PEOLE in tabular form. The time 
system is UTC. It is hoped that these elements will prove useful to other in- 
vestigators. 
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Table 4 

Mean Orbital Elements for PEOLE in 1971 



TIME (MJQ) 

A (METERS) 

E 

1 (DEGREES) 

W (DEGREES) 

N (DEGREES) 

M (DEGREES) 

40966.462777718 

7C09597. 65794 

O 1 63180 72 

15.005041 

284 • SO 234 5 

329.983986 

272.467983 

40969. 59 85 87 93 t 

7009496,24795 

.0 1 6431391 

15.003364 

326.512804 

?0«. 1 36B96 

72.322104 

40960. 535451 350 

7009065.18692 

.0 16760 6 87 

15.00248 0 

109.920721 

231 .911210 

73.006024 

4098 1. 527499936 

7008989.375 1 5 

.3 1 6763530 

1 5.004165 

122.806151 

224.997559 

323 .522267 

40906. 527499936 

7008760 .49327 

.0 I 6439371 

15.005973 

108.881498 

19C .138038 

347. 85991 3 

40996.691 770B07 

70CS434. 12950 

•0 1 6328953 

1 5.0068B5 

324.96781 l 

1 1Q.2570P0 

1 94 . 275*50 

41 00 7. 73 9 095 800 

7000 153.84 322 

.3 J 6731696 

15.005039 

109,975350 

42.207647 

0? .77421 B 

41010. 02 3229 1 33 

7000093. 9S8 1 3 

• D 16596132 

15.005453 

150. 325923 

2f , 707355 

326* 8206^ 5 

4101 3. 875349697 

7009C29. 01813 

.0 l 6407595 

15. 006312 

190. 634667 

359.409315 

46 . 245626 

41014. 720844675 

7008023.29546 

.0 1 6364Q75 

15. 005303 

201.625101 

353.52184S 

235.94925? 

4101 7.61 3 t 36546 

7007979.37785 

.3 1 62331 29 

15. 004977 

240.499837 

333.352564 

1 02. 61C 42B 

41021. 56*91 B947 

70Q7931 .57C96 

'• 0 l 6213049 

1S.003C86 

293. 677443 

305. 802025 

20*165382 

41027.302175809 

7007828.91640 

.316499448 

15.005339 

10.04561 1 

2 65.786622 

Z? . 22! |06 

41 031 .51 CSC9222 

7C07731 .82794 

.0 1 6680768 

1 5.003418 

65. 226442 

236.432323 

15C .427399 

41033. 74 93 0S556 

7007704.26232 

.31 6682622 

15.003252 

94.608071 

220.803376 

2 1 2. 8 32C62 

4k 030.665972223 

7007579.79792 

.0 1 6520497 

1 5.003250 

158.651671 

I 86 .502305 

1 66 .7 31 429 

4t C4 1 . 66 5972223 

7007510.57854 

.0 t 6315559 

13.004321 

198. 360467 

165.572664 

323. 0*361 1 7 

4 1 04 t . 6659 72 223 

7007511 .13861 

.31 6316708 

15.003988 

198.348222 

165.572256 

323.906145 

41 044. 07388RA58 

7007430.93191 

.0 16207604 

1 5.006617 

233 . 469847 

148.782859 

208 .59196,2 

41044.073888858 

7007424.3624 1 

.0 16191831 

15.007851 

233.531730 

14R. 771502 

208.532473 

41 055. 70 7638889 

7C07C 12.27887 

.3 16473809 

1 5.007148 

26. 056669 

67.617400 

?4«. 053476 

41060. 555995330 

7006790.47464 

.0 1 661 1 171 

15.005030 

89.41 6683 

33.785076 

297.275156 

41061. S67870335 

7006748.1 1580 

.316600013 

IS, 004057 

102. 51 7254 

26.718476 

296 • 2 C 2 66 3 

41 067, 566643483 

7C06553. 07094 

.0 16315249 

15.005819 

161.265120 

344.851391 

262.722309 

4107C.S05451 365 

7006f 29.290 1 5 

• 0 161 478R9 

15.006345 

223 • 426074 

324 .341752 

102.501 169 

41076.937997651 

7C 06457.96380 

.016109947 

15.007621 

306. 899043 

279.467015 

222.535129 

41080.521330984 

7C06401 .82527 

.01 6309 65C 

15.008408 

354.695 80 9 ! 

254.4591 90 

2M . 580 107 

41083. 737141 183 

7006364.99685 

.3 t 6489700 

15.005768 

37. 046092 

232.0 10002 

140.419192 

4| 067. 71 26C41 32 

7006300 .71825 

.016589230 

15.006053 

00. 935690 

204.245469 

113. 253548 

41093. 540972223 

7006207.76799 

.0 16338391 

15.004546 

1 65. 19884? 

163.546331 

253.613233 

41 09 8. 70 76389Q9 

70061 66.71 B70 

.0 1 6050260 

15.0C4907 

234 . 02 6188 

127.486426 

102.151093 

41 105. 624305556 

7006075.61479 

.016123877 

15.006421 

325.94439? 

79.204999 

288 . 09069 1 

41 106. 595466001 

7C 0 60 56 - 4 0 5 25 

.0 16174233 

15.0 06 81 3 

339.906843 

72.427259 

70.272310 

41113. 49 5305556 

7005929.8673 1 

.01 6526Q77 

1 5.007457 

70. 74241 2 

24.221344 

IQ?, 680570 

41 11 3.495305556 

7005929.71787 

.01 65C0672 

15.004844 

70.683663 

24.220347 

190. 734C0 1 

41 119.624305556 

700S83I .39543 

.3 1 6397 1 85 

1 5,005168 

150.857124 

341 .457069 

I 16.358894 

41 124.624305556 

7005772 .62825 

.0 l 6046345 

1 5. 002939 

217. 156717 

306.558090 

158.1 14630 

41128. 74 93 0555 6 

7005713.31352 

.01 5960 1 38 

1 5.003575 

272.469959 

277.770430 

210.203634 

41 134. 665972223 

7005619. 615B0 

,0 16200890 

J5.00543Q 

351.795978 

236.450749 

103.661178 

41 134,665972223 

70056 19. 47976 

.0 1 6191 483 

IS. 005 103 

351 • 826308 

236. 461012 

103.626994 

41 14 3.666574022 

70055 42.42048 

.01 6502691 

15.002993 

109. 867984 

173.586867 

257.670506 

41 147.583240689 

70055IT.31B55 

•0 1630B187 

15.001 108 

161 ♦ 163059 

146.232117 

?RCi .272045 

41 151.54 3310151 

7005491 .635 1 3 

• 0 1 6C6681 1 

16.001788 

213* 816710 

1 10.576584 

173.BZ4C12 

41 206.456724507 

70048 14.716 1 2 

.0 1 60 0 7 8 96 

15.007528 

221 • 425675 

95.064661 

210.917568 

41232. 15902780 l 

7004522. 1 5697 

.016007106 

15*004148 

201 . 475885 

275. 513842 

241 . 771 183 

41237. 796035231 

7004429.07921 

.01 5879B33 

15.003299 

277. 554388 

236. 1 25368 

90 . 0 9 1 C 6 1 

41243. 354455969 

7004350.33496 

.01 6135797 

14.988362 

351.952758 

197. 311994 

240. 10451 4 

41248.98 9837927 

7004243.731 43 

.0 1 6164420 

14.981067 

65. 708509 

1 57.930460 

03 , 1 97075 



CONCLUSIONS 


The new center of mass coordinates presented in this paper in combination with 
our previous (Marsh, et al., 1973) results give strong ties covering the Eurasian 
continent and northern and southern Africa. The mean orbital elements of 
PE OLE. determined fro m a combinat ion of laser and Minitrack data, have al- 
ready proven useful in geodynamic investigations (King Hele, 1973; Wagner, et al., 
1974). All of these results demonstrate the value of properly coordinated co- 
operative programs. 
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